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Abstract 

We propose two dynamical models with delay taking advantage of 
their complex dynamics for information processing tasks. The first 
model incorporates coupled delayed dynamics of multiple bits, which 
is shown to have desirable properties as an encryption scheme. The 
second model is a single binary element with delayed stochastic tran- 
sition, which presents a resonance behavior between noise and delay. 



1 INTRODUCTION 

Complex behaviors due to time delays are found in many natural and 
artificial systems. Some examples are delays in bio-physiological controls 
|2[ , and signal transmission delays in large-scale networked or distributed 
information systems (See e.g. ||, |]]). Research on systems or models with 
delay has also been carried out in the fields of mathematics |5|, g|, artificial 
neural networks @, |j, and in physics H [R], 11]. This series of research 



has revealed that time delay can introduce surprisingly complex behaviors 
to otherwise simple systems, because of which delay has been considered an 
obstacle from the point of view of information processing. 

In this paper, however, we actually take advantage of this complexity 
with delayed dynamics and propose two models. The first model incorporate 
delayed dynamics into a new model of encryption. The second model shows 



1 To appear in the proceedings of 1998 International Symposium on Nonlinear Theory 
and its Applications (Nolta'98) September, 1998, Crans-Montana, Switzerland 
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a resonance behavior between noise and delay. Through these two models, 
we explore the possibility of applications for delayed dynamics. 

2 THE ENCRYPTION MODEL 

The encryption process is identified with a coupling dynamics with various 
time delays between different bits in the original data. We show that the 
model produces a complex behavior with the characteristics needed for an 
encryption scheme. 

Let us now describe the encryption model in more detail. S(0) is the 
original data of N binary bits, whose ith element Sj(0) takes values +1 or 
— 1. The delayed dynamics for the encryption can be specified by a key which 
consists of the following three parts: (1) a permutation P generated from 
(1,2,3, -,N), (2) a delay parameter vector r which consists of N positive 
integers, and (3) number of iterations of the dynamics T. Given the key 
K = (P, t, T), the dynamics is defined as 

Si{t) = (-l)XS pi {t- Ti ), (1) 

where pi and Tj are ith element of P and r, respectively. (If t — T{ < 0, we 
set t — Ti = 0.) In Figure 1, this dynamics is shown schematically. The state 
of the ith element of S(t) is given by flipping the state of the pjth element of 
S(i — Tj). Thus this dynamics causes interaction between N bits of the data 
in both space and time. The encoded state S(T) is obtained by applying 
this operation of equation (1) iteratively T times starting from S(0). 




t - Xj t - Tj t 



Figure 1: Schematic view of the model dynamics. The state of the zth 
element of S(t) is given by flipping the state of the pith element of S(t — Tj). 

We investigate numerically the nature of the delayed dynamics from the 
perspective of measuring the strength as an encryption scheme. First, we 
examine how the state S(t) evolves with time. In Figure 2 (A), we have 
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shown an example of encoded states with different T using the same P and 
r for a case of N = 81. To be more quantitative, we compute the following 
quantity as a measure of difference between two encoded states at different 
times t and if. 



A typical example is shown in Figure 2(B). We note that the dynamics 
of our model has a property of occasionally very similar, but not exactly the 
same, states appearing indicated by sharp peaks in the figure. Except around 
these particular points, however, we generally obtain rather uncorrelated 
encoded states (i.e., Y ~ 0) with different iteration times. This is a desirable 
property of the model as an encryption scheme: the same state can be 
encoded into uncorrelated states by changing T. 

Next, we investigated the effect of a minor change of P and r on the 
model dynamics. Starting with the same initial condition, we evaluate how 
two states S(t) and S'(t) are encoded with slightly different P and P', re- 
spectively, by computing 



A representative result is shown in Figure 3(A). The same evaluation with 
r and t' is shown in Figure 3(B). These graphs indicate that if we take 
sufficiently large T, the same state can evolve into rather uncorrelated states 
with only a slight change of P and t. This again is a favorable property in 
the light of encryption. It makes iterative and gradual guessing of P and r 
in terms of their parts and elements very difficult: a nearly correct guess of 
the values of P and r does not help in decoding. 

With these properties of the model, an exhaustive search appears to 
be the only method for guessing the key. Even if one knows N and T m ax, 
the largest element in r, one is still required to search for the correct key 
from among (Nl)(T max ) N combinations and to guess T. The commonly used 
DES (Data Encryption Standard) employs 2 56 bit keys ]12[ . We can obtain 
a similar order of difficulty with rather small values of N and T max ; for 
instance, N ~ 11 and T max « 10 PL 

There are different methods possible for using this model for a secure 
communication between two persons who share the key. One example is 
that the sender sends a series of encoded data in sequence for the interval 




i=0 
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(3) 
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Figure 2: (A) Examples of encoding with the model dynamics from an 
initial state to T = 500 and T = 1000 with P and r. (B) Example of the 
correlation Y between encoded states at different time steps evaluated by 
equation (2). Cases with tf = 1000 are plotted with the initial state and P 
and t the same as in Figure 2 (A). 



between T and T + r ma x (or longer). The receiver can recover the original 
data from this set of encoded data by applying a reverse dynamics with the 
key. In a situation where the data sent is a choice out of multiple data sets 
known to the receiver, the receiver can run the encryption dynamics to the 
entire sets with the key for case matching. 

3 RESONANCE BETWEEN NOISE AND DE- 
LAY 

The second model is a stochastic binary element with delayed transition 
which is schematically shown in Figure 4. The state of the element X(t) 
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P' = 



4 33 62 49 47 26 20 7 21 
65 72 57 73 67 75 25 34 40 
60 12 70 10 14 43 38 17 29 
51 46 77 81 52 45 31 54 39 
76 71 30 44 28 16 13 8 56 

1 78 22 19 24 18 61 74 59 
66 64 53 27 80 63 68 5 
[55] 23 11 9 79 58 35 42 15 
50 41 48 37 36 6 2 32 69 
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Figure 3: Example of the correlation X evaluated by equation (3) between 
two states encoded by slightly different (A) P and P', and (B) r and r 1 . 
The difference of P' from P, and r' from r are indicated by boxes. The 
initial state and P and r are the same as in Figure 2. 



at time step t can take either —1 or 1 and the transition probabilities are 
indicated by the arrows. We further include delay in the model, which makes 
the formal definition of the model as follows: 

P(l,t + 1) = v X(t-r) = -l, 

= 1-q X(t-r) = l, (4) 

P(-l,t+l) = q X(t-r) = l, 

= 1-P X(t-r) = -l, (5) 

where P(s,t) is a probability that X(t) = s. Hence, the transition probabil- 
ity of this model depends on its state at r steps earlier. In this sense, this 
model is a special case of " delayed random walks" [|ll], [14| except it can take 
only two states. 

Let us now discuss our motivation for considering such a model. It has 
been shown that a simple stochastic binary element can show a resonant 
behavior with a external oscillating signal with an appropriate choice of 
transition probability. This phenomena is termed as "Stochastic Resonanc" 
(see e.g. (l^]) now widely investigated. In our model here, we have replaced 
an external oscillatory signal with delay and its associated oscillatory be- 
havior. Hence, we can expect some form of resonance can be seen with the 
model between the oscillatory dynamics due to noise and that due to delay. 
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A preliminary simulation study shows this is indeed the case. We have 
fixed q = 1 — q = 0.5 and r = 10, and varied p. The simulations are 
done up to 20000 steps. In Figure 5, where we have plotted examples of 
residence time histograms in the — 1 state and dynamics of X(t), we can see 
a resonance phenomena as a height of the peak at the period equal to the 
delay. If we tune the noise (p) appropriately, we observe the height of the 
residence time histogram peak reaching the maximum value. This is shown 
in Figure 6 [16]. 

With this property, this model can be used to stochastically encode in- 
formation. We have built a preliminary application to encode static pictures 
using this model, which will be reported elsewhere. 



4 Summary 

We have proposed two models which include delay. These two models indi- 
cate that complex delayed dynamics can be used for information processing 
tasks rather than necessarily being obstacle. More thorough investigation of 
each model is presented elsewhere, and concrete applications of these models 
are currently being devloped. 




Figure 4: Schematic view of the stochastic model dynamics. 
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Figure 5: The residence time histogram and dynamics of X(t) as we change 
p. The values of p are (A) p = 0.005, (B) p = 0.04, (C) p = 0.2. 
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Figure 6: A plot of peak height with varying p 
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